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This work has been stimulated by the quenching of spontaneous emission in an optically thin high- 
density plasma reported by Suckewer and coworkers [Phys. Rev. Lett. 60,1122(1988)]. Our work is 
based on the collision-induced coherence of two decay channels along two optical transitions. The 
quantum interference of pumping processes creates the dark state and the more atoms are pumped 
in this collision-induced dark state the stronger the suppression of the spontaneous emission. The 
efficiency of this suppression is quantified by putting it in comparison with the spontaneous emission 
on the ultraviolet transition which proceeds in a regular fashion. The branching ratio of these 
two(visible and ultraviolet) transitions is introduced as the effective measure of the degree of the 
suppression of the spontaneous emission on the visible transition. Our preliminary calculations show 
that a significant decrease of the branching ratio with increase of electron densities is reproduced in 
the theory. More experiments and more elaborated theories are needed to conclude about the role 
of quantum coherence and interference in such manifestly incoherent media as plasmas. 



I. INTRODUCTION 



An excited free atom decays into a vacuum emitting spa- 
ially isotropic radiation in the spectrum of frequencies which 
jas Lorentzian shape with a bandwidth proportional to the 
Einstein A coefficient [jpj. The decay rate is a combination 
fundamental atomic constants and strongly resistant to 
>a modification. However, such modifications are possible, 
syfor instance when the atom is no longer free and placed in 
^a high-Q wavelength-size cavity. The decay of the strongly 
^confined atom has been the subject of intensive theoretical 
--.find experimental research^. Here, the rate is enhanced 
Sii)r suppressed through the modification of the density of 
^electromagnetic modes in the neighborhood of the resonant 
,_}requency[3i- 5:]. Using a three- level atom driven by coherent 
^fields is another way to control the spontaneous emission [6| 
. J^nd observe narrowing of spectral linewidth compared to the 
^<4iatural linewidth 0, H| . Four-level driven atomic configura- 
tions offer yet another level of control allowing, for example, 
. .Spectral line elimination and full spontaneous emission can- 
celation HOI]. 

The topic of above-mentioned studies has been focused on 
individual transitions in atoms and addressed the modifica- 
tion of one or another rate of spontaneous radiative decay. A 
number of studies consider more involved schemes character- 
ized by correlated emission simultaneously from two atomic 
upper levels. Such type of interference of decay channels has 
been exploited, for instance, by Harris in his proposal of las- 
ing without inversion in fTlj ] . In contrast to most studies, we 
consider not the direct interference of decay channels, but 
rather the interference of incoherent pumping rates (namely, 
collision-induced interference). In its turn, this interference 
induces the correlation in spontaneous emission. We propose 
using branching ratio R as the measure of this correlation. 
The term branching ratio here stands for the ratio of two 
spontaneous decay rates from two upper levels |a) and |6) to 
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FIG. 1. Five- level configuration of working levels. The three-level 
subsystem \a) — \c) — \b) (inside the oval) is the core for the study 
of the suppression of the spontaneous emission on the visible tran- 
sition. Here 7 a and 7& are spontaneous decay rates and r a and rb 
are collision-assisted pumping rates. For simplicity they are set 
pairwise equal: j a = j b = j vis , r a = r b = r vi3 . The \a) o \d) is 
the ultraviolet transition along which the spontaneous emission 
proceeds in a regular fashion with rate yuv- The pumping rate 
r e to state \e) with simultaneous equal decay rates 7 e to states 
I a) and \b) is important as these processes pump the collision- 
induced dark state, as explained in the text. The rate serves 
to close the system. Note that pumping rates along all dipole- 
allowed transitions are present in plasmas. The particular value 
of the rate depends on the energy of the transition according to 
the Boltzman distribution. 



two lower levels |c) and \d) in a five-level atom of Fig[TJ 
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where for definition we formulate the branching ratio in 
terms of decay rates on visible, j V i s , and ultraviolet, juv> 
transitions. Here, we refer to particular transitions in order 
to keep a close link to the previous experiments on measure- 
ments of changes of branching ratios. These experiments 
were performed in plasmas in Princeton by Suckewer and 
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coworkers [121 ] . They serve as the main motivation and ref- 
erence point of our study, although our results are applicable 
to more general schemes. 

The branching ratio is closely related to total spectral- 
line intensities (in photons) for the two transitions, I V i S and 
luv ■ This relation is simply the consequence of the fact that 
the intensities are proportional to corresponding decay rates 
(for schemes considered so far): I V i S oc 7„j s and Ijjv oc luv- 
In fact, this proportionality allows us to consider the inten- 
sity as alternative operational definition of the decay rate 
which differs from the pure decay rate by a factor. The ac- 
cess to this factor can be, however, a difficult experimental 
task. This potential experimental problem disappears when 
we consider the branching ratio. In those cases where the ra- 
diation is emitted from same level (s) the factors are the same 
for both transitions and therefore cancels out in the ratio in 
Eq. ((J). So, the operational definition of the branching ratio 
can be introduced as 



R = 



LUV 



(2) 



Moreover, for some cases, including ours, the operational 
definition in Eq. © is very natural and probably the only 
possible definition. The problem is that we consider simul- 
taneous decay along two (visible) transitions. This decay is 
characterized by more than one decay rate. This problem 
does not appear for the ultraviolet radiation, the decay of 
which happens along one transition and is therefore charac- 
terized by a single decay rate. However, we cannot formulate 
the branching ratio in the form of Eq. ([I]); and instead, we 
characterize the spontaneous emission process by the opera- 
tional definition given by Eq. @. 

For optically thin samples and decays from common up- 
per level, the equality between two definitions Eqs. flTJ and 
([2]), i.e. the equality between the ratio of line intensities 
and the ratio of the corresponding decay rates, points to the 
way of measurements of decay rates. Thus, from measure- 
ments of the intensity ratio of two lines, the ratio R and 
one of the spontaneous decay rates can be deduced knowing 
the other. Such spectroscopic technique is rather popular in 
plasma physics. Particularly attractive is the fact that the 
branching ratio is as immune to environmental conditions as 
the decay rate is. For instance, the ratio does not depend on 
the electron density and temperature. In the experiments of 
it is reported that for some transitions the branching 
ratio becomes the function of density (concentration). They 
interpreted this observation as quenching of spontaneous- 
emission coefficients for visible transitions. Our study is 
devoted to formulation of the problem in quantum-optical 
terms, and our model is based on the five-level scheme (see 
Fig. [TJ, and the effect is caused by collision- induced quan- 
tum coherence between two upper levels. 

The best way to formulate our main result is to compare 
two cases. The first one corresponds to no coherence be- 
tween upper levels a and b. Then, the branching ratio for 
visible and ultraviolet transition found from the operational 
definition $2$ reads 



as found in the following sections. Here the difference in 
the factor of two from original definition (fTJ) arises due to 
doubling the visible decay rate (here i a — lb = Ivis) by tak- 
ing into account decays from both upper levels (with equal 
rates). In contrast, the rate on the ultraviolet transition is 
not doubled because the spontaneous decay originates from 
single level a. Apart from this factor, the equation ^ stands 
in line with regular expectations. In particular, this result 
is unable to explain the concentration dependence of R re- 
ported in (l2T |. 

The second case realizes when the (maximal) coherence 
between levels a and b is induced by electron collisions. 
Then, as we shall show in the following sections, the branch- 
ing ratio under steady-state conditions is given by 
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where, for simplicity, the result is evaluated in the limit of 



large pumping rates: 



, r e > 7„ ls , juv ; and fa is the ap- 
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propriate cross-sections of electronic excitation. Since all in- 
volved pumping rates depend on the electron concentration, 
the branching ratio also becomes concentration-dependent. 
The higher the concentration the faster the pumping rates 
and therefore the lower the branching ratio. This relation im- 
plies quenching of spontaneous emission on the visible tran- 
sition. 

In the following sections we shall show that the effect of 
quenching is explained by creation of the collision-induced 
(i.e. induced by pumping rate r V i S ) dark state as a linear 
combination of upper states a and b. This dark state does 
not emit radiation and the more atoms we pump into this 
state the better suppression of spontaneous emission we ob- 
serve. Indeed, as shown below, by pumping more atoms 
in the dark state via increasing r e with concentration we 
achieve more efficient quenching. In contrast, the emission 
on the ultraviolet transition does not depend on the coher- 
ence between levels a and b and is linearly proportional to 
population of a level only. By formulating the branching 
ratio R, all common dependencies cancel out and only the 
purified asymmetry of responses on the visible and ultravi- 
olet transitions is left. Here, the branching ratio shows up 
as a valuable measure of the coherence-induced quenching of 
spontaneous emission on the visible transition. 

Overall, in this chapter we report the scheme with sponta- 
neous emission cancellation resulted from incoherent pump- 
ing of atoms into the collision-induced dark state. As quan- 
titative measure of the degree of the cancellation we propose 
to use the operational definition of the branching ratio given 
in Eq. ©. Our scheme provides a possible physical expla- 
nation for the quenching effects reported in [12j. 



II. BASIC SET OF EQUATIONS 

We study the set of equations of motion governing tempo- 
ral evolution of the density matrix elements for the five-level 
scheme shown in Fig. [TJ The key ingredient of this scheme is 
the three-level V -type system consisting of two upper states 
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| a), | b) connected to the lower state \c) by dipole allowed 
transitions. We first consider this three-level subsystem and 
then add two states \d) and |e). We assume that two transi- 
tions | a) — |c) and \b) — |c) have close transition frequencies 
and therefore upper states decay to same continuum of vac- 
uum modes with decay constants j a and 7&. These constants 
as well as details of the decay process follow from the Hamil- 
tonian 

k 

+ hY / 9k b} e l( ^~" k)t \b)(c\a k (5) 

k 

+ H.c. 

describing interaction between the atom and the reservoir of 
vacuum oscillators, each of frequency Vk {k here represents 
both the momentum and polarization of the vacuum mode) . 

Here are the coupling constants between the fc-th vac- 

uum mode and the atomic transitions from \a) and |6) to |c). 
a>k{a\) is the annihilation (creation) operator of a photon in 
the fc-th vacuum mode, which obeys conventional bosonic 
commutation rule [afc,ajj = 8kk' ■ 

Electron collisions transfer atoms (also ions as in the 
plasma experiments in Ref. [Hj]) between states in the man- 
ner very much like that in the process of interaction with an 
optical field. As a consequence, the Hamiltonian appears in 
a similar form, as 

Vr = PaMa){c\ + p bc n\b)(c\ + H.c. (6) 

In contrast to coherent optical fields, individual collisional 
events are not correlated with each other and therefore the 
effective field f2 is incoherent. Note that spectrum of this 
incoherent field is very broad so that it covers both upper 
levels simultaneously, and therefore one and the same field 
drives both transitions. These transitions are characterized 
by dipole matrix elements p ac and pb c and transition fre- 
quencies Lo ac and ujbc- The incoherent field has <5-like corre- 
lations at different times, i.e., 

(Q*(t)Q(t'))=RS(t-t') (7) 

These correlations are not expected to cover the entire range 
of frequencies from —00 to +00. It is sufficient that they 
are at least approximately valid in the vicinity of both res- 
onances and cover the frequency separation of two upper 
levels. 

The total interaction picture Hamiltonian for the three- 
level subsystem is the sum of two terms introduced above, 

H = V 7 + V R (8) 

Derivation of equations of motion for the density matrix ele- 
ments from the total Hamiltonian is a straightforward task. 
It is based on the scheme well developed in quantum op- 
tics, see for instance [Hj], where main ingredients are the 
Wigner-Weisskopf approximation and the generalized reser- 
voir theory More details, particularly on the 6 system 
under consideration, can be found in our recent paper [151 ]. 



With this reference to prior works we skip the derivation and 
jump directly to the equations of motion. They read 

Paa = ~{r a + "fa)Paa + r a Pcc ~ \^P \fr^Xb{Pba + Pab) (9) 

pbb = ~(n + lb) Pbb + npec - \p\fr a rb{pba + Pab) (10) 



Pec = {r a + la)Paa + (n + -fb)Pbb - (r Q + n)p C c (11) 

+ PV r an{pba + Pab) (12) 

Pab = ~\{ r a +r b +la+ lb)Pab ~ i&Pab (13) 

+ ^PV r a r b(^Pcc ~ Paa ~ Pbb) (14) 

Pea = ~^( 2r a + r b+ la) Pea ~ ^PV r o.nPcb (15) 

Pcb = -\{^r b + r a + Jb)Pcb ~ \p\/r a nPca (16) 



Here r a .b = 2{p 2 acbc /h 2 )R are pumping rates of atoms to 
upper states \a) and \b) induced by collisions with electrons, 
see Fig. [T] Detuning A = uj ac — u>i> c is supposed to be small 
as compared to the optical frequency and should be of the 
order of pumping rates r a or or less in order to make 
quantum interference effects observable. 

Terms containing products of pumping rates appear due 
to the interference of two optical transitions which is induced 
by collisions with electrons. These terms are central to our 
discussion. In order to emphasize the role of this collision- 
induced interference, we assume that decay channels do not 
interfere. The interference terms are accompanied by the 
p factor. This factor is the normalized scalar product of 
corresponding dipole matrix elements: 



\pacPbc\ 

According to its definition, the alignment factor takes the 
value equal to 1 for parallel dipole moments, —1 for antipar- 
allel, and for orthogonal. Intermediate values on [—1,1] 
segment are also possible. Maximal coherence corresponds 
to parallel or antiparallel dipole moments, while zero coher- 
ence corresponds to orthogonal dipole moments. These two 
extremes of maximal and minimal coherence deserve special 
attention. 

It is not difficult to proceed with modeling the system 
shown in Fig. [T] and supplement the three-level subsystem 
with two additional levels. For further considerations we 
shall need only equations for populations and one equation 
for the polarization on \a) — \d) transition. The upper state 
|e) is subject to decay with rate 7 e and pumping from state 
|c) by electronic collisions with rate r e . So, for the population 
we immediately obtain 

Pee = r e (Pcc ~ Pee) ~ 2-lePee (18) 

Of course, we could derive this equation from first princi- 
ples by writing the Hamiltonian of interaction of the atom 
with the continuum of quantum oscillators to model the 
decay and interaction with the incoherent electron field to 
model the pumping process. Since no interference effects are 
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expected on this transition, this derivation is no more than 
a trivial exercise. The intuitive appearance of this equation 
stays in contrast to the detailed consideration of the V -type 
subsystem where the nontrivial appearance of interference 
terms requires cautious analysis. 

Similar equation can be written for the lowest state of the 
system, \d). Here the decay with rate juv from \a) state 
pumps the \d) state. The electron pumping is modeled by 
rate rjjv- So, in accordance with the picture in Fig. [T]we 
write 



pdd = ruv(p cc - Pdd) + lUVPa 



(19) 



We now complete the set of equations by writing the equa- 
tion of motion for the polarization on \a) — \d) transition: 



Pad = —-^{lUV + ruv)Pad 



(20) 



We can also write down equations of motion for the other 
polarizations involving states |e) — \d) . However, we do not do 
this because these equations are not used in further analysis. 
For more general considerations, if necessary, they can be 
composed by following directly the picture in FigfT] Note 
that no interference terms appear. 

To complete the model we should account for decays and 
pumpings between the ^-subsystem and \d) and |e) states. 
The subset of equations (|9"))- (fTrJ]) are modified in the follow- 
ing way. There is one incoming term for the population of 
|6) state: +j e Pee- One incoming and one outcoming term for 
Paa ■ +lePee and -"fuvPaa- For p cc we get two additional 
pumping terms, +r e (p ee - p cc ) and +r uv (p dd - p cc ). The 
equation of motion for the two-photon coherence p a b should 
be supplemented with decay term ^uvPab- One-photon po- 
larization p ac also decays faster due to additional decay term 
\{r e + Tvv + luv)Pac- Finally, in the equation for pb c we 
should also write additional decay term \{r e + rjjv)pbc- In 
the next section we summarize all these changes and after 
implementing a few simplified assumptions write down mod- 
ified equations. 



III. DRESSED-STATE ANALYSIS 

In this section we consider properties of the V-type subsys- 
tem \a) — \b) — |c) in the dressed-state picture. Before doing 
this we first prefer to make some simplifications. Assume 
equal dipole moments on \a) — |c) and \b) — |c) 8 transitions: 
Pac = Pbc- This assumption immediately yields the equality 
of decay constants: j a — jb- We also get equal pumping 
rates, r a — r . Also, let the two transitions have equal fre- 
quencies, so that two-photon detuning A is zero. 

The upper level |e) is auxilary with respect to the rest 
of the system. It serves to model the external pumping 
of two upper levels of the V -subsystem. In plasmas this 
pumping is due to thermal redistribution of population from 
all levels which are higher with respect to \a) and \b). In 
order to accurately account for details of this process we 
would have to consider a particular level scheme for a given 
atom/ion. However, our immediate goal is not the one-to- 
one description of a particular experiment but rather the 



proof-of-principle theoretical study of relevant coherent ef- 
fects. Guided by simplicity we eliminate level |e) from our 
five-level scheme. For that purpose, we assume fast decay 
rate 7 e with respect to pumping rate r e . Then, from the 
equation (|18p we can write for the population p ee in steady- 
state: 



Pee ~ — Pc 

le 



(21) 



With these simplifications in effect and additional terms 
from the discussion in the end of the previous section, equa- 
tions (P|)- (IT^1) for the V -subsystem become 

Paa = ~(r v is + Ivis + lUv)Paa + (r v is + ~^)Pcc ~ V r visPab 

(22) 

- pr V isPab (23) 

- r e + r uv )p cc (24) 



Pbb = -{r V is + ivis)pbb + (r V is + y)Pc 

Pcc ij^vis 4" 'Yvis^)(,Paa 4~ Pbb} (^Tvis 
+ ^Pr vis Pab + WvPdd 

^iuv 1 

Pab = "{rvis 4- Ivis 4 ^~)Pab + -jV r vts(^Pcc - Paa - Pbb 



(25) 

5bfc) 

(26) 



with new notations: j a = 75 = ^ V is and r a = Tb = r V i S - Here 
we write down only equations relevant for the study in this 
section. Equations for one-photon polarizations p ca and p c b 
will be considered in the next section. 

Coherent effects in the V-subsystem appear due to the 
interference of pumping channels along \a) — |c) and \b) — |c) 
transitions. Actually, it is the same (incoherent) pumping 
process that drives both transitions. The presence of the 
coherence shows up in most natural way in the dressed-state 
basis. Here, the pumping field dresses atomic states \a) and 
I b) yielding the collision-induced dark state: 



\D) 



r b\a) ~ Vr^\b) 
\/r a + r h 



and its orthogonal partner - bright state: 



\B) 



r a \a) + y/n\b) 

V r a + n 



(27) 



(28) 



Here we formulate the dark and bright states for general 
relationship between pumping rates to emphasize that our 
treatment is valid beyond the simplification assumption of 
equal dipole moments. When r a — rb the expressions become 



\D) = -j=(\a) \b)) 



\B) = ^=(\a) + \b)) 



(29) 



(30) 



In the following we shall see that the bright state has direct 
relationship to the rate of correlated spontaneous emission 
from bare upper states |a) and |6). 
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Since Schrodinger description is impossible for the given 
scheme of decay and pumping rates, the formulation of equa- 
tions of motion in terms of the dressed variables requires in- 
troduction of corresponding density matrix elements. They 



1 . 

PDD = -^{fiaa + Pbb ~ 2pab) 



PBB = -j{Paa + Pbb + ZPab) 



(31) 



(32) 



This understanding comes from the analysis of the equation 
(j3U). We shall see that more atoms in the dark state makes 
the emission weaker. The first two terms in equation (134)) 
stand for the depletion of the dark state and are therefore 
undesirable for our purposes. [Note that steady-state polar- 
ization pbb contributes as little as ("fuv /r V i S )pDD-} This 
depletion is counterbalanced by the gain associated with the 
last term. The large value of this term requires fast pumping 
rate that is guaranteed by inequality ([37| . 

With inequality (I37[) we get approximate solutions of equa- 
tions They are 



PDB ~ -^{Paa — Pbb) 



(33) 



where we assume definitions of the dark and bright states 
given in equations ([2"§f and ([3"Uf and real two-photon polar- 
ization p^. The last is valid assumption for zero two-photon 
detuning. If necessary, it is rather straightforward to derive 
the density matrix elements for general relationship between 
pumping rates and for nonzero two-photon detuning. 

Given equations of motion in bare-states basis, (|2"2"j) . (|23l) 
and (f2"6")l . and dressed-states definitions (|3"Tj) - (|3"3"f for the den- 
sity matrix elements, we reformulate these equations as 



PDD 



-IvisPDD 



luv 



(pdd + Pdb) + r e p c 



(34) 



Pbb 



-(2r„ 



)PBB- 



luv 



(pBB+PDB) + (2r vis +r e )p cc 
(35) 



PDB = -(r v is+lvis + 



-^)pdb-^1uv(pdd+Pbb) (36) 



Our ultimate goal is the comparison of intensities of spon- 
taneous emissions on the visible transition on one hand and 
on the ultraviolet transition on the other. As we shall derive 
in the next section, these intensities are linearly proportional 
to the population of \B) state and |o) state, correspondingly. 
Both these populations can be obtained from the solution of 
equations (j34l) - (j3"6"|) and inverted conversion formulas (|3Tj) - 
([53"| . It is not difficult to obtain exact steady-state solutions 
for all three dressed-states density matrix elements in above 
equations. However, we prefer to make yet another simplifi- 
cation in order to bring physics in the clearest possible form. 

We assume that pumping rates are much faster than de- 
cay processes, particularly faster than the fastest decay rate 
7(/v- Note that the spontaneous emission rate juv is two 
orders of magnitude greater than decay rate "f V i S on the vis- 
ible transition. Clearly, this difference is attributed to the 
cubic dependence of spontaneous decay rates on the transi- 
tion frequency. In its turn, the high efficiency of pumping 
processes arises from frequent collisions of ions with free elec- 
trons (note that the electron concentration N e exceeds the 
value of 10 18 cm -3 ). Finally, we formulate our assumption 
in the form of strong inequality 



(37) 



This inequality is the key to understanding the effect of 
suppression of spontaneous emission on the visible transition. 



Pdd 



Pbb ~ (1 



luv 



4r„ 



-)Pc 



Pdb ~ Pc 



(38) 



(39) 



(40) 



where we also used juv 3> Ivis- All these solutions are 
expressed in terms of p cc which remains unknown quantity 
until we solve all equations of motion for the closed five- 
level system. However, the explicit knowledge of p cc is not 
necessary for our study. 

Solution (|3"5|) shows that the dark state is populated only 
due to the pumping via auxiliary upper state |e). By com- 
paring solutions given by formulas (|38j) and (|39j) and using 
inequality (f3"T|). we conclude that pdd Pbb, so that the 
population of upper states \a) and \b) is dominantly concen- 
trated in the dark state. This asymmetry means that most 
atoms are trapped in the non-emitting state \D). As we shall 
see shortly, the spontaneous emission on the visible transi- 
tion is proportional to pbb\ and therefore basing on the just 
derived strong inequality pdd 3> 7bb, we conclude that 
the emission is strongly suppressed. This suppression is to 
be put in comparison with the regular (unsuppressed) emis- 
sion on the ultraviolet transition, which is regulated by the 
amount of population in the \a) state. Overall, the degree of 
suppression is quantitatively estimated by the corresponding 
branching ratio, as the ratio of two emissions, see Eq. ([2]) 
and derivations in the next section. 

The population of the \a) state as well as p b and p a b can 
be found by inverting definitions (|31|) -(|33 |) . So, in terms of 
bare states we get 



1 r e 
t:(pbb + Pdd + 2p DB ) ~ tt^i 



Pbb = t;(pbb + Pdd - 2 pdb) ~ n 

2 Zj uv 



1 T e 

Pab = t:(pbb - Pdd) ~ Pc 

2 2j uv 



Pcc 



(41) 



(42) 



(43) 



This state is the state of maximal coherence in the sense that 
\Pab\ ~ \fp~aapb~b- This coherence is of collisional nature and 
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its high degree is the reflection of the high efficiency of the 
pumping mechanism (r vis ). So, we deal here with almost 
pure states. 

In this section we performed the analysis of the V - 
subsystem in terms of dressed states and introduced the 
important approximation (|37|) which is the key to under- 
standing of the effect of suppression of spontaneous emission 
on the visible transition. In the next section we proceed 
with derivation of the emission rates on visible and ultravi- 
olet transitions and put them in comparison in terms of the 
branching ratio defined via Eq. @. 

IV. SPONTANEOUS EMISSION SPECTRA 

The idea is to calculate the spectrum S V i S (uj) of spon- 
taneously emitted photons on the visible transition, i.e. 
through | a) — |c) and |6) — \c) channels. Then, we simi- 
larly calculate the ultraviolet spectrum Suv(w) of photons 
emitted in the \a) — \d) channel. Integrating these two spec- 
tra over the frequencies in the vicinity of the corresponding 
emission peaks and normalizing to the energy of one photon, 
we thus obtain two spectral-line intensities 

I vis = (huyis)" 1 J dujS vls (uj) (44) 

luv — (f^Juvy 1 J dwS V i s (u>) (45) 

Taking the ratio of these intensities we formulate the branch- 
ing ratio defined in Eq. ([2]) , finally in terms of relevant pump- 
ing rates and decay constants. 

Spontaneous emission spectrum of a two-level atom be- 
came already a textbook subject in quantum optics, see for 
instance [13]. Here the characteristic three-peak spectrum 
emerges in an elegant way as a result of application of the 
quantum regression theorem. The latter allows us to ex- 
press the two-time correlation function in terms of one-time 
correlation functions and thereby significantly simplifies the 
calculation procedure. 

Spontaneous emission spectra of coherently driven three 
and more level atoms were the subject of interest in numer- 
ous papers, see for instance [16j . Again, the quantum regres- 
sion theorem is of great help. However, even in the case of 
three levels the calculation procedure quickly becomes cum- 
bersome and requires numerical analysis on the final stage. 
Details can be found in the original paper (l6| . In our case 
the system is driven by an incoherent field that makes the 
situation a little simpler. Furthermore, applying the so far 
discussed approximations and simplifications, we can make 
the problem even analytically tractable with transparent fi- 
nal result. We even bypass the direct application of the 
quantum regression theorem. 

First, we take the definition of the spectrum (valid for sta- 
tionary processes) as the Fourier transform of the autocor- 
relation function of the second order of the scattered electric 
field: 

S(oj) = -Re dTe- luJT (E-(r,t) ■ E + (r,t + r)) (46) 

Jo 



This autocorrelation function is simply related to the 
normally-ordered product of polarizations taken at instants 
of time separated by a positive time delay r. This relation- 
ship follows from the well-known formula 

E+(r,t) = ^nx(nxd)pW(t-r/c) (47) 

47T£oc 2 r 

connecting the electric field E^ + ) with polarization P( + ). 
Here, n is the unit vector in the direction of observation, d 
is the unit vector along the atomic dipole moment, r is the 
observation point, measured from the position of the atom, 
and loq is the polarization frequency. More details can be 
found, for instance in Ref. |l6[. The resultant formulas are 

S vis (u) = —^Re / dre-^{p- is {t) ■ P+ S (t + r)) (48) 
k Jo 

fV,4 poo 

S uv (u:) = ^^Re dre^ {P^ v (t)-P+ v {t+r)) (49) 
* Jo 

Here C is an unimportant constant, same for both tran- 
sitions. P^J and P^J (P^j/and Pf)\)) are negative and 
positive frequency parts of the polarization induced on the 
visible(ultraviolet) transition. According to our scheme they 
are defined as 

Pvisi*) = PacOac{t) + Pbc&bc(t) (50) 
Pvisi* +T) = pca<Tca{t + t) + p cb (T cb (t + t) (51) 

Similarly we define the negative- and positive-frequency 
parts of the polarization on the ultraviolet transition: 

PUV® = Pad<?ad{t) (52) 
P UV ( t + T ) = PdaCfda (t + T) (53) 

In the following we assume p oc = p ca = p& c = p c & and 

pad = pda- 

Operators of atomic transitions <Jij are defined in usual 
way: aij = They are particularly relevant for the 

calculation of multi-time correlation functions. Note that 
these quantities are quantum-mechanical operators. Their 
quantum-mechanical average has simple relation to the ele- 
ments of the density matrix, namely 

(a ij (t))=Tr[U + (t)a i j(0)U(t)p(p)] 

= Tr[a lJ (0)U(t)p(0)U+(t)} (54) 
= Tr[a ij (0)p(t)} 

Since we need to calculate two-time correlation functions, 
it is instructive to write down equations of motion for rele- 
vant operators and then find solutions to the initial-value 
problem. The form of these equations coincide with the 
equations of motion for respective density matrix elements. 
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This can be checked by direct derivation from the Heisen- 
berg equation of motion with the Hamiltonian Thus we 
get 



1 



1 



r (3r vis + r uv +r e +7„ is +7 £ /v)(Tac - 7; r vis<?bc (55) 



CTfcc = -~(3r vis + r uv +r e + j vis )a bc - \^r ms a ac (56) 



cfad = --^{r V i s + r uv +r e + j vis + Juv^ad (57) 

where we used the assumption r a = r b = r V i S . Note that 
(Tji = ajj and the equations of motion for the conjugate 
operators arise by taking the Hermite conjugate of the right- 
hand sides of Eqs. (|55l) - (|57|) . The above equations form the 
closed system and allow simple solutions. These solutions 
acquire even simpler form when we apply our key inequality 
([571) . Thus, for the quantities of interest we get 



c C a(i 



,(t)+<7«*(t)] (58) 



(59) 



where r = (2r vis + r uv +r e ). 

Everything is ready for the formulation of the desired re- 
sult. After substituting solution (|5"5]l in the expression (|5 1 [) 
and correspondingly (|59[) in (|53[) . the two-time correlation 
functions of the polarizations read 

(Pvis(t>Pvis(t+r)) = p 2 ac e- roT ((a aa (t)+a bb (t)+a ab (t)+a ba (t)) 

(60) 

(Puv(t) ■ P$v(t + t)) = p 2 ad e-^ +r ^(a aa (t)) (61) 

The Fourier transform of the above expressions converts 
the exponential decays in time to the Lorentzian spectra in 
frequency: 



C^li S p 2 ac r 

7T UJ 2 + 7'n 



2p BB {t) 



(62) 



S uv (w) 



C ^uvP 2 ad r uv +r v 



+ {ruv + r vis y 



;paa(t) (63) 



Here we take only the real part of the integral and use the 
rule (|54[) for calculating quantum-mechanical averages, and 
also the definition of the population of the bright state (|3"2"j) . 
Note also the equality p ab = p ba valid for the zero value of 
the two-photon detuning, which is the case considered here. 

Note here that it may seem that the thus obtained spectra 
can, in principle, become a function of time, copying the time 
dependence of populations pbb and p aa . This dependence is 
a signature of a nonstationary process. Therefore the spec- 
trum of these solutions is to be defined not by Eq. (|46[) but 
with the more general formula which is derived and discussed 
for instance in Ref. [13j. However, in our case the interest 



is in stationary emission corresponding to steady-state solu- 
tions for pbb and p aa . Therefore, Pbb arid p aa become time 
independent. In Fourier space they are simply two station- 
ary Lorentzian peaks. Spectral-line intensities follow as the 
frequency integral over these peaks: 



I vls - h 1 Cujl is pl c 2p B B 



IUV = ft 1 CuJ UV pl d Paa 



(64) 



(65) 



The branching ratio is the ratio of two above quantities, 
as regulated by definition ([2]): 



R = 



Jvis %PB 



B 



IUV Paa 



(66) 



where we used the definition of the spontaneous decay con- 
stant 



1 4w„ J 



4"7T£o 3frc 3 



(67) 



and similar formula for the ultraviolet transition, and can- 
celed out common prefactors. 

This ratio is the main formula of our chapter. It re- 
flects the role of coherent effects in the spontaneous emis- 
sion along the visible transition. A small population in the 
bright state on the background of relatively strongly pop- 
ulated | a) state implies the suppression of the spontaneous 
decay. Such asymmetry becomes possible only when the co- 
herence p ab between the two upper states \a) and \b) acquires 
an appreciable value. Here, the coherence is induced by colli- 
sions and its large value is guaranteed by the inequality (34) , 
which means that decay processes destroying the coherence 
are of little importance. Quantitatively, the degree of the 
suppression is related to the value of R in absence of the 
interference terms which becomes simply the double ratio of 
decay constants, as we show below. 



V. SUPPRESSION OF THE SPONTANEOUS 
DECAY 

The branching ratio given by Eq. (po]) is now the subject 
of the analysis with respect to the experimental observations 
in Ref. [12|]- First of all we use steady-state values of pop- 
ulations given by solutions and (|4ip in Eq. (|6"6"]l and 
obtain 



R = ivi 



l 



The common factor p cc has been canceled out. According 
to the experimental observations in Ref. (T^ |. the branching 
ratio was measured as the function of electron density. One 
order of magnitude increase in the electron density resulted 
in one order of magnitude decrease in the branching ratio. 
Let us see how such dependence emerges from the above 
equation. 

The decay constant 7^ depends nontrivially on the fre- 
quency and the dipole moment, see Eq. (|67p . Within the 
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range of change of the experimental conditions, no depen- 
dence of the transition frequency (i.e., no line shift) on elec- 
tron density was found. Some theoretical efforts were made 
in Ref. [17| to take into account the effect of the screening 
of the atomic potential by the surrounding plasma. This ef- 
fect can, in principle, modify the Coulomb field experienced 
by the valence electron and, thus, significantly change the 
transition probability (i.e. dipole moment). However, this 
change becomes of an appreciable value only for concentra- 
tions much higher than that used in the experiment. So, we 
conclude that no change takes place in the value of j V i S when 
concentration varies. 

The observed changes of the branching ratio R in Ref. 
[l2| can be attributed mainly to the dependence of colli- 
sional rates r e and r vis (and, of course rjjv) on the density 
of free electrons. Apparently the higher the density the more 
frequently the collision events occur and therefore the larger 
the rates. More precisely, we assume linear dependence of 
the collisional rates on the concentration (r, oc N) and write 
for a constant electron temperature 

n (N) = k,N (69) 

where i = e, vis, UV, and ki is the proportionality coefficient. 
Typical range of concentrations where the effect occurs cov- 
ers the region 10 18 — 10 19 cm~ 3 . Applying dependences given 
by Eq. (|69[) to Eq. (|68[) . we perform simulation and present 
obtained results in Fig. [5J We thus illustrate main result of 
our study - the sensitive dependence of the branching ratio 
on the concentration. The range of the change is of the same 
order as was observed in the experiments. 

For completeness we analyze the dependence of the colli- 
sional rates (and therefore, the branching ratio) on the elec- 
tron temperature of plasma. This dependence, even when 
present, cannot be deduced from the experiments in Ref. [12[ 
where the electron temperature was estimated as constant 
under operating conditions and equal to « 5eV. However, 
in the set of related experiments on measurements of the 
branching ratio in high-density plasmas of CIII performed 
in a later work in Ref. [l8j . the temperature varied consid- 
erably and, therefore, we can expect a well pronounced de- 
pendence of the branching ratio on the temperature. Thus 
for the concentration 0.7 x 10 18 cm -3 , the temperature was 
5.7 eV; while for higher concentration 2.6 x 10 18 cm -3 , the 
temperature increased to 9.3eU. 

It is instructive to analyze the simultaneous electron den- 
sity and temperature dependences. Free electrons in plasmas 
obey the classical Boltzmann distribution thus the number of 
particles dN (in a unit volume) within the range of energies 
E + dE yields 

dN = N^ JJ ± m e-^VEdE (70) 
So, for collisional rates we can write 



— zkl -E- 
ri (N,T) = 2NkiJ—e-** (71) 



where ki is a cross-section. Then the coefficient ki is given 

by 

— IlkT E - 

ki=2 H^ e '^ (72) 

The branching ratio was measured in Ref. [18] for five 
concentration/temperature points. 

VI. SIMULATION RESULTS 

We have performed numerical simulation to demonstrate 
the dependency of branch ratio on electron density. For pur- 
poses of plasma diagnostics, we simulate the situation with 
or without the coherent effects due to electron impacts, de- 
pending on the relative orientation of the dipole moments of 
optical and UV transitions. 

Here we assume that j V i 8 = 1, so that all the other 
parameters are normalized by "f V i s . The collision- induced 
incoherent pumping rates are dependent on electron den- 
sity; here we use rjjv = 0.001 x N e , r e — 0.1 x N e , and 
r v = 0.3 x N e . We considered three different UV decay 
rates with yuy = 0.1,1,5 to investigate the dependence of 
branching ratio on electron density. 

Fig[2](a) corresponds to small UV decay rate juv = 
0.1-f vis . Fig[2J(b,c) correspond to larger UV decay rate 
luv = 1)5 respectively. In each figure, the top red curve 
is for p = 1, the bottom blue curve is for p = — 1, and the 
dashed line is for p = 0. 

Through these nine cures, we notice that, at large electron 
density, the dependence of the branching ratio on density for 
different p factors behave in a similar way and trend to 0.5. 
In the special case with p — 0, this process gets to limit much 
more quickly. It is also clear that with larger UV decay rate, 
the dependence of the branching ratio on density for different 
p factors also behaves very similarly, but at a lower decay 
rate. At both sides of the extremely low or high electron 
density, there clearly limits. Outside this range, the electron 
density has no effect on the branching ratio. 

Obvious conclusion from the comparison of curves is that 
the increase of temperature greatly suppresses variations in 
the branching ratio with concentration. The relatively small 
slope of the curve could be the reason or one of the rea- 
sons that yielded the very weak concentration dependence 
and led the experimenters in Ref. [18| to the conclusion of 
the absence of the dependence of the branching ratio on the 
concentration. 

We thus complete the comparison with the experiments. 
On our way to the main result expressed by formula (I68[) 
we made two important assumptions. One is the condition 
of maximal coherence. The other is the domination of col- 
lisional rates over all relevant decays, see Eq. (|3T|) . It is 
instructive now to see what happens if we relax the first as- 
sumption and consider the opposite case - the case of no 
coherence. Simply set p factor to zero; that is equivalent to 
setting pah = (at least in the steady-state). Given inequal- 
ity (|37[) . two upper states \a) and |6) are populated equally 
in steady state. In the case of no coherence, the population 
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Log(N ) 



FIG. 2. The dependence of the branching ratio on logarithm of 
electron density, (a) For j uv = 0.17^, (b) For j uv = 7«<»,(c) For 
7uu = 57i;i s . Three lines in each figure correspond to p — 1,0,-1. 



of the bright state is simply one half of the sum of p aa and 
p bb . So, we get 



Paa ~ Pbb ~ (1 + tt)Pcc 
2r 

With these solutions in Eq. (|6"6")l. we finally obtain 



R 



■ 7« 



luv 



(73) 



(74) 



for the branching ratio for the visible and ultraviolet transi- 
tions in the case of no coherence between two upper states |o) 
and 1 6). This value is considerably larger than the branching 
ratio obtained for the case of maximal coherence, see Eq. 
(jnHJt- Moreover, this value does not depend on collisional 
rates and therefore on concentration. The comparison of 
branching ratios calculated for these two cases demonstrates 
the key role of coherence in interpreting experimental results 
obtained in Ref. (T^ |. 

It can be noticed that the branching ratio Limit, at low 
and high electron densities is 



1 + 2-^- N e ^0 

f? — ) luv c 

luv 



(75) 



With no coherence excited by electron collisions, the ratio 
is given by 



R 



jJs/uv + 2{jv + r v )) 
luvilv +r v ) 



(76) 



VII. DISCUSSION 

So far we discussed the suppression of spontaneous decay 
on the visible transition. This suppression is due to the two- 
photon coherence induced by the interference of collisional 
processes along two overlapping optical transitions \a) — \c) 
and 1 6 — \c). In order to quantify the degree of the sup- 
pression we introduced the branching ratio as the ratio of 
total spectral-line intensities for the two transitions. Here, 
one transition is the visible transition (actually, simultane- 
ously two transitions) of interest and the other transition is 
the reference ultraviolet transition. This branching ratio cal- 
culated for the case of maximal two-photon coherence, see 
Eq. (jB"5)) . is compared to the branching ratio evaluated for 
the case of no coherence, see Eq. (|74)) . The degree of the 
coherence-induced suppression can be deduced as the ratio 
of these two branching ratios, i.e. 



, 4 1 

R = luv [ — I 



(77) 



The formula is valid as long as inequality (|3T[) holds. The 
faster the collision rates the stronger the interference of the 
optical channels and the larger the degree of the suppression. 
Physically, the suppression of the spontaneous decay on the 
visible transition is associated with putting most atoms in 
the non-emitting dark state, while only a small fraction of 
the atoms left in the bright state participates in the emission, 
as comparison of populations given by formulas (35) and (36) 
illustrates. 

The important message of this study is the demonstra- 
tion of the possibility of coherent effects in plasmas, where 
the coherence is induced by the interference of incoherent 
processes. Here the collisions of free electrons with ions rep- 
resent these incoherent processes, and they indeed happen 
rather frequently for concentrations as high as 10 18 cm" 3 . 
When they dominate over relevant decay rates, the two- 
photon coherence becomes of substantial value that leads 
to efficient suppression of the spontaneous emission on the 
visible transition. This suppression was registered experi- 
mentally and reported by Chung, Lemaire, and Suckewer in 
Ref. [l2( . Simultaneously, no suppression was registered for 
the spontaneous radiation on the ultra-violet transition. So, 
the effect shows up when these two emissions are compared 
in terms of the branching ratio. With our theoretical model 
we can reveal and explain these features. Moreover, we are 
able to explain sensitive dependence of the degree of suppres- 
sion on concentration of free electrons. It is interesting that 



10 



the dependence on the concentration only appears due to the 
collision-induced coherence and disappears in the absence of 
the coherence, as comparing Eqs. (|68p and (|74p shows. 

Motivated by the comparison in our study with the ex- 
periment, we draw attention mainly to the role of quantum 
coherence in suppression of the spontaneous decay. An in- 
teresting question is whether this coherence can cause the 
opposite effect, namely the enhancement of the spontaneous 
emission. The answer is positive, provided we modify a little 
our model. Let us eliminate the pumping rate r e via auxil- 
iary upper state |e), while keeping other pumping and decay 
rates intact. Effectively, the five- level model reduces now to 
four levels. Equations (l34l) - (|36|) have simple solutions. They 
become even simpler under our main assumption that the 
pumping rate r V i S dominates over the other rates. The pop- 
ulation of the dark state decays to zero while the bright state 
is equally populated as the lower state |c): 

Pbb ~ Pec (78) 

Since the dark-bright coherence is as small as 
— {luv / ^r V i S )pBB, the population of |o) state is given 
by 

11 

Paa ~ -^PBB ~ ^Pcc (79) 

Taking the ratio as in Eq. (|66j) . we formulate the branch- 
ing ratio for the visible and ultraviolet transitions: 

R = 4^- (80) 

luv 

We get twice as strong emission as in the case with no co- 
herence, compared to Eq. (|74[) . This doubling is the largest 
that we can get as the enhancement factor for the rate of 
spontaneous emission. 

Of course, this case with zero value of r e is not realiz- 
able in plasmas where all pumping rates scale with energy 



according to the Boltzman distribution (|T1 [) and cannot be 
controlled from outside. However, this case demonstrates 
another facet of the quantum coherence. Here, as in the 
case of suppression, the coherence is again of large (maxi- 
mal possible) value, but now with opposite sign: p a t, w \p C c- 
Physically, all populations of upper states are concentrated 
in the emitting bright state \B), while the non-emitting dark 
state is empty. Therefore, the emission becomes as strong as 
it can be. 

In conclusion, we suggested the model for demonstrating 
quantum interference effects in multi-level systems and fol- 
lowing [l2| proposed the branching ratio defined according 
to formula ^ as the measure of suppression/enhancement 
of spontaneous emission in multi-level configurations. We 
applied the scheme to the interpretation of experimental re- 
sults reported in Ref. [l2j and studied in detail the effect of 
suppression of spontaneous emission on the visible transition 
as opposed to the absence of the suppression effect on the 
ultraviolet transition. 

In particular, we revealed the characteristic dependence 
of the branching ratio on the concentration of free electrons: 
an order of magnitude increase in the concentration corre- 
sponded to one order of magnitude decrease in the branch- 
ing ratio. To our knowledge, this interpretation is the first 
attempt in explaining experimental results in high-density 
plasmas from the viewpoint of quantum interference. More 
experiments and more elaborated theories are needed to con- 
clude on the role of quantum coherence and interference in 
such manifestly incoherent media as plasmas. 
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